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Abstract. Coherence rephasing is an essential step in quantum storage protocols 
that use echo-based strategies. We present a thorough analysis on how two adiabatic 
rapid passages (ARP) are able to rephase atomic coherences in an inhomogeneously 
broadened ensemble. We consider both the cases of optical and spin coherences, 
rephased by optical or radio-frequency (rf) ARPs, respectively. We show how a 
rephasing sequence consisting of two ARPs in a double-echo scheme is equivalent to 
the identity operator (any state can be recovered), as long as certain conditions are 
fulfilled. Our mathematical treatment of the ARPs leads to a very simple geometrical 
interpretation within the Bloch sphere that permits a visual comprehension of the 
rephasing process. We also identify the conditions that ensure the rephasing, finding 
that the phase of the optical or rf ARP fields plays a key role in the capability of the 
sequence to preserve the phase of the superposition state. This settles a difl:erence 
between optical and rf ARPs, since field phase control is not readily guaranteed in the 
former case. We also provide a quantitative comparison between 7r-pulse and ARP 
rephasing efficiencies, showing the superiority of the latter. We experimentally verify 
the conclusions of our analysis through rf ARP rephasing sequencies performed on the 
rare-earth ion-doped crystal Tm'^"'':YAG, of interest in quantum memories. 
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1. Introduction 

Very recently, adiabatic rapid passages (ARP) liave drawn tlie attention of the optical 
quantum storage community [I112], as they offer the possibility of rephasing atomic 
coherences. In inhomogeneously broadened materials for quantum storage purposes, 
the fact that the atomic coherences dephase is the key property that allows the storage. 
At the same time, bringing the coherences back into phase is a necessity for retrieving 
the stored information (see, e.g. [3]). The coherences to be rephased can be both of spin 
or optical nature. Even though optical ones are needed to interact with the information- 
carrying photons, conversion to spin coherences is preferred in between the capture and 
the retrieval stages of the storage protocol in order to take profit of their longer lifetime. 
In such schemes, coherences of both types dephase and need to be rephased by optical 
and radio-frequency (rf) means [Ij. This is the case of echo-based memories. Protocols 
based on electromagnetically induced transparency carefully eliminate the stage where 
the information is stored in the optical coherences and allow direct storage in the spin 
coherences [5l[6] by means of optical control pulses. In this case, only rf rephasing is 
needed j2]. However, this technique is not best-suited for inhomogeneously broadened 
media due to its narrow-band characteristics. Neither matter the specific characteristics 
of the protocol nor the nature of the coherences, typically, these are rephased by the 
application of a vr pulse. However, disadvantages concerning vr pulses have been pointed 
out. One is the need of high powers for the optical or rf fields, which scales as the 
protocol bandwidth squared. Another one, concerning optical coherences only, is the 
unavoidable distortion of the tt pulse as it propagates through the sample because of light 
absorption [7]. This makes that the pulse area is no longer vr after some threshold depth 
in the sample. The third disadvantage, related to the previous one in the optical case, 
is the high sensitivity of the technique to spatial intensity variations. These drawbacks 
have made researchers explore ARPs as an alternative. 

An ARP consists of a pulse whose frequency is chirped through a range that 
typically goes from much lower to much higher frequencies than the ones in the 
inhomogeneous width one wishes to rephase. The amplitude of the field can either 
be varied as well or kept constant (see figures [3](a) and[3]^b)). 

Adiabatic rapid passages have been extensively used in nuclear magnetic resonance 
(NMR) through decades. Initially, in the early 1960 's, it proved to be a means for 
intensity- insensitive and frequency-selective adiabatic spin inversion [8]. Transfer to 
the optical domain came some years later in the experiments by Loy [9l[T0]. A rich 
investigation on pulse shape optimization (amplitude and phase modulation), from 
analytically [HHIH] to numerically [T5[fT8] - [2Tj derived or proposed methods, opened 
the way to increase the intensity insensitivity and/or the frequency selectivity of the 
pulses. 

In 1987, Kunz approached the question of using ARPs for spin rephasing [22]. He 
highlighted the fact that rephasing is not possible with a regular ARP. Put in words 
more suitable for the present context of quantum information storage, the reason he 
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stated is the following: the passage of the ARP, leaves the spins with a phase that is a 
function of the spin transition frequency of each particular atom. This phase needs to be 
compensated for in order to enable spin refocusing. From Kunz [22] and other research 
groups (see [23],|2l] and references in [25]) raised a variety of more complex pulses, 
generally combining regular ARPs or half ARPs as building blocks, that provided phase 
compensation. However, rather early, Conolly and coworkers realized that just adding 
a second identical ARP is enough to compensate for the phase induced by the first, 
and they proved it experimentally [26]. Here, words need to be carefully chosen not 
to mislead the reader. The phase compensation occurs in the reference frame that 
rotates at the instantaneous frequency of the rephasing pulse field, and still further 
clarification is needed for the term phase compensation to be exact since that frame 
changes abruptly from the first to the second ARP. In the laboratory frame, the term 
"phase compensation" is not appropriate. The picture is more complicated and we 
will discuss it in detail in the section that follows. In any case, coherences can be 
rephased indeed by the application of two ARPs. The pulses need not be consecutive, 
which allows for storage time. Rf realizations in the context of NMR can be found in 
refs. [T8 |[271[28] . More recently, interest in ARPs has risen in the quantum information 
community. Specific implementations for coherence rephasing have been demonstrated 
both in the rf [Tl|2l[29] and optical [30] domains in the context of broadband quantum 
memories. Nevertheless, recent applications of ARPs do not limit to their rephasing 
capabilities. Their use for optical-to-spin coherence transformation and back has been 
theoretically studied in ref. [3l], showing that they preserve the collectivity of the 
superposition state of an atomic ensemble (i.e., the phases of each of the components 
involved in a Dicke state). Other kinds of adiabatic pulses, which satisfy the adiabatic 
condition (to be defined later) but which do not involve frequency chirps have also 
been developed. The Stark-chirped rapid adiabatic passage (SCRAP) compensates the 
difficulty of chirping short (nanosecond) optical pulses by inducing a time- varying shift in 
the atomic transition through a far off-resonance adiabatic optical pulse. Another time- 
delayed close-to-resonance adiabatic pulse performs the population transfer between two 
optical states [32l|33]. As another example, the stimulated Raman adiabatic passage 
(STIRAP) relies on the adiabatic variation of the amplitude of two detuned and time- 
delayed pulses to generate a population transfer between two spin states by optical 
means (see [31] and references therein). 

Two articles provide enlightening analysis of how rephasing by two ARPs happens. 
These are refs. [25], by Hwang and Shaka, and [35], by de Graaf and Nicolay. The 
latter authors use the most spread out approach in NMR, that is, 3D geometrical 
representations of the trajectory of the effective magnetic field (the control vector, in 
a language better suited to the present article) and that of the magnetization vector 
(the Bloch vector) under the effect of the former. Hwang and Shaka proposed a much 
more concise approach by choosing a matrix treatment of the effect of each ARP. 
Then, a rephasing or any other sequence is constructed as the product of the matrices 
corresponding to the building blocks. This way, conclusions can be drawn very easily. 
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However, Hwang an Shaka limited their analysis to one particular type of ARP-based 
sequences, of interest in NMR (intended to eliminate the deleterious effects of the huge 
water resonance in the quality of NMR spectra). 

In this article, we propose a matrix treatment of the rephasing sequence involving 
two ARPs, in the Bloch sphere formalism. We will revisit and extend the analysis 
of the matrix computation for one ARP only, already undertaken in a previous work 
by some of the authors of the present one [1]. Later, we will make use of the ARP 
matrix to build up the matrix associated to a rephasing sequence. We will see that this 
matrix approach allows a simple understanding of how the rephasing works and why 
two ARPs are needed to rephase an inhomogeneously broadened ensemble of atomic 
coherences. A geometrical interpretation derived from the matrix treatment provides a 
visual comprehension of the process. We will draw comparisons with vr pulses and regular 
TT-pulse-based echo sequences in terms of rephasing efficiencies. We will also analyse the 
role of the rephasing field phase, as well as the phase preservation characteristics of 
the ARP-based rephasing sequence, particularly relevant for quantum memories. At 
this point, a significant distinction between rf and optical rephasing will be highlighted, 
originated in the different capabilities each technology provides for controlling the (rf or 
optical) field phase. We will present experimental verification of the critical conditions 
that enable rephasing. Radio-frequency spin rephasing experiments are performed on a 
rare-earth ion-doped crystal, a thulium-doped YAG, a material that has been actively 
studied as a candidate for quantum memories [36H11]. 

The paper is organised as follows: In a theoretical section we will deal with the 
derivation of the matrix expression for an ARP from the Bloch equations, the matrix 
construction of a two-ARP rephasing sequence, the comparison with vr pulses and 
considerations on phase preservation. The experiments are presented in the section 
coming after. Finally, some concluding remarks will be given. 

2. Theory 

The ARP matrix will be derived thanks to successive frame changes. The first one is 
well-known and uses the rotating frame within the rotating wave approximation (RWA). 
It allows an accurate and consistent relation between the excitation field phase and the 
atomic state phase. The second frame change takes the tipping control vector as a 
polar axis. It provides a simpler interpretation of the different angles in our geometric 
interpretation. 

2.1. One adiahatic rapid passage 

2.1.1. The ARP matrix In this section we aim at obtaining a matrix expression for an 
ARP. We will disregard decoherence effects. In what follows, we will use the Bloch sphere 
representation for the dynamics of a two-level system, whose states will be denoted \a) 
and 1 6), in interaction with an external oscillating field A{r,t). The field can be either 
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electric or magnetic, depending on the nature of the ARP, optical or rf, respectively. 
In such representation, the two-level system state is expressed by the Bloch vector B 
defined as 

^ Pab + Pha \ 



\B] 



\ 



KPba 
Pbb - 



-pab) 
Paa 



(1) 



where pij represents the element {i\p\j) of the density matrix operator p and the notation 
[ ]r indicates that the vector coordinates are given in the laboratory reference frame R 
of cartesian axes u, v, w. We express the field as 



A(r, to + t) = Ao(t) cos {coot + 0(t)) 



(2) 



where Ao(t) is a real vector, uq and (p{t) are real quantities and to is the central instant 
of the ARP. In the Bloch sphere representation, the parameters that define the field 
compose the control vector in the following way 

/ 2n{t) cos {uot + (f){t)) \ 




3(to + 1) 



J R 



\ 



(3) 



/ 



where Q{t) is the Rabi angular frequency of the field and Uab = {Eb — Ea)/h is the 
angular frequency of the transition between the states \a) and \b) of energies Ea and 
Eb, respectively. In the case of an optical field, we have Q{t) = /x^^ ■ Ao(t)/^, with 
H^jj = {a\er\b) the electric dipole moment of the transition (— e is the electron charge). 
We have chosen the relative phase of the states \a) and \b) so that is real. Thus, 
Q is real as well. In the case of a magnetic field, we have Q = 7-Bi, where 7 is the 
gyromagnetic factor and Bi is the component of the oscillatory field perpendicular to 
the quantization axis. 

In the Bloch sphere representation, Schrodinger's equation takes the simple form 



B{t) = ~m^B{t), 



(4) 



in the absence of decoherence. Equation (jl]) describes the precession of B around f3. 
The last statement might mislead the reader to think that the dynamics of B are rather 
simple. As a matter of fact, f3(t) is a vector that oscillates at frequency uq + (pit). 
Therefore, a vector precessing around it describes a complicated trajectory. 

The picture becomes simpler if we change to a reference frame where the unit vector 
f3 = /3/\l3\ is static. We will perform the transformation in two steps. First, we will 
change from the laboratory frame R to the reference frame R' that rotates around the 
vertical axis at frequency uo + (j){t) with phase 0(0) at t = 0. R' is known as the rotating 
frame. The matrix 

^ cos{uot + (pit)) sin(co'o^ + 0(^)) 
Ciit) = -sin(wo^ + 0(t)) cos{ujot + (pit)) 






1/ 



(5) 
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performs the frame change from R to R' through the operation [B{tQ + = 
Ci{t)[B{tQ + This change of reference frame is equivalent to applying the 

transformation pabito + 1) H- p'ai,{tQ + t) = Pabito + 1) exp{-i{uot + 0(t))), pu H- p'^ = pu 
(with i = a,b) to the density matrix. As regards the control vector, its coordinates in 
the frame R' are given by 

f n{t) [1 + cos{2{ujot + ^{t)))] \ 
^{to + t)]= -n{t) sm{2{uJot + (j){t))) , (6) 

y UJab J 

Taking the RWA, the above expression reduces to 

/ n{t) \ 

kto + t)] = , (7) 



R' 



\ 





^ab 



J 



The Bloch vector dynamics in frame R' are ruled by 



R' 



B{t) - (wo + ^{t))w X B{t) 

- H 

X B{t) - {uo + <j)it))w X B{t) 
I3{t) X B{t) 



(8) 



with a new control vector f3 defined as 



\ 



\ 



(9) 





A-<p{t) J 

and A = Uab — Wq. The dynamics of B in frame R' given by ([8]) are much simpler 
than in frame R. Indeed, in R' and under the RWA, the fast oscillatory behaviour of 
the control vector has been ruled out. What is left is just its smooth variation as Q{t) 
and (f){t) evolve during the ARP. The smoothness of that variation is ensured by the 
adiabatic condition to be described shortly. The Bloch vector now precesses around the 
slowly varying axis f3(t) with angular frequency 

In an ARP, the frequency of the external field is varied through a wide range, from 
values much smaller than Uab to values much larger than Uab (or inversely, depending 
on the sign of 0). In an ideal positively chirped ARP, the instantaneous detuning of the 
external field from the atomic frequency, A — <p{t), varies from oo to — oo. At the same 
time, Q(t) stays bounded, i.e. \fl{t)\ ^ oo. Therefore, in the course of the ARP, the 
unit vector /3{to + 1) somehow goes from w' to —w', with w' the vertical axis of frame 
R' (w' = w). The exact trajectory and instantaneous angular velocity of [/3(to + t)]R' 
will depend on the specific way Q{t) and (f){t) are varied during the ARP, as well as 
on the detuning A. The simplest case, occurring for a constant Q, a linearly chirped 
frequency (0(t) = rt) and zero detuning halfway through the ARP (A = 0) is depicted 
in figure [H The precession of the Bloch vector is also represented. 
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4>i-T/2) 2 



Figure 1. (Colour online) Bloch sphere scheme of an ARP satisfying A = and 
|17(±r/2)/<^(±T/2)| < 1. We have chosen (j){t) = rt with r > and n{t) = n. 
At t = — T/2, the control vector (3 points upwards. The control vector slowly starts 
turning upside down. During the lapse -17/(2(/«(-T/2))T < t < VL / {2^{-T / 2))T , 
the control vector stays off the vertical axis. This is a very small fraction of T given 
the assumption mentioned above. Halfway through the passage, the control vector 
is in the equator. At the end of the passage, the control vector points downwards. 
All throughout the passage, the Bloch vector precesses around the control vector, 
describing the cones shown in the figure. On each snapshot, the Bloch vector is 
represented at an arbitrary position within the cone. 



We mentioned above that a second change of reference frame is useful to further 
simphfy the trajectory of the Bloch vector. In this case, we will change from R' to a 
reference frame R" where the vertical axis w" is instantaneously parallel to /3(to + 
This can be achieved through the matrix 



/ cos^(t) -sin^(t) \ 



C2{t) 



1 
y sine(t) cosOit) J 

which ensures [B{tQ + t)]/?" = C2{t)[B{tQ + t)]/?'. The angle 6{t) is defined as 

A - . VL{t) 



cos^(t) 



1/2 



sin6'(t) 



fi(t)2 + (A-0(t))2 

and is sketched in figure [U In R" , f3(to + t) takes the form 

/ \ 



1/2- 







V 



n{tf + (A - 

Thus is static. The dynamics now become 



1/2 



(10) 



(12) 



R" 



B{t) - e{t)v X B{t) 



J R' 



(5{t) - e{t)v' X B{t) 



(13) 
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The dynamics simplify under the adiabatic approximation. It assumes that B essentially 
precesses around /3, i.e. 



9{tf < (3{tf 
or, equivalently, 



(14) 



n{tY + (A - 



-1 3/2 



< 1 



(15) 



The adiabatic approximation assumes that the precession rate (generalized Rabi 
frequency) is much faster than the control vector tipping rate. The condition is usually 
found in literature expressed as 

— > 1, 16 

r 

obtained for the worst-case scenario (A = 0) occuring for a chirp such that f2 = and 

0(t) = r (linear chirp). 

If f|T5|) is valid, the dynamics reduce to 



m 

Integration yields 



R" 



I3{t) X B{t) 



(17) 



[Bito + h)]R. = U{h, Q [B{to + ta)]R 

' COS x{tb,ta) - sin xih, to) \ 

U{tb,ta) = Smxitb,ta) COS X{tb,ta) 

^0 1 



n{ty + (A - <j){t')y 



1/2 



dt'. 



(18) 
(19) 

(20) 



where x(4, to) represents the total precession angle 

We have now gathered all the necessary elements to build a matrix expression for 
an ARP. The matrix associated to an ARP connecting the initial and the final states is 



M, 



ARP 



Ci(T/2)-iC2(T/2)-i[/(T/2, -T/2)C2(-T/2)Ci(-T/2), (21) 



with T the duration of the ARP. Due to the inclusions of Ci and C2, Marp acts on the 
reference frame R. It is worth noting that Ci and C2 are only necessary to be known 
at the instants of the beginning and end of the ARP. 



Securing coherence rephasing with a pair of adiahatic rapid passages 



9 



2.1.2. Far off-resonance initial and final conditions As we have said above, an ideal 
ARP is the one where /3 departs from w' at tQ — T/2 and reaches —w' at tQ + T/2, in the 
case of a positive chirp (from lower to higher frequencies). We can simply understand 
that by analysing the case of an ARP-driven population inversion experiment. The 
Bloch vector is initially parallel to ±w'. The ARP will most efficiently drive the Bloch 
vector all the way to the =pti>' direction only if f3 is initially parallel (or anti-parallel) to 
B and remains so all throughout the passage. In practical realizations, this condition 
is most closely satisfied if 

n{±T/2) 

A - 0(±T/2) 

The detuning should be larger than the Rabi frequency or, in other words, the initial 
and final excitations should be far off-resonance. The expression has to be satisfied for 
any atomic detuning A within the inhomogeneous broadening Pinh- 



< 1. 



(22) 



If (122!) is satisfied, the matrices C2(— T/2) and C2(T/2) reduce to 



Co 



T 
2" 



Co 



( 



\ 










1 





(23) 



-1/ 



Therefore, the mathematical treatment of the 
significantly. 

In what follows we will assume ( l22l) is valid. 



ARP matrix fl2T]) will simplify 



2.1.3. Eigensystem In section [2.1.11 we have deduced an analytical expression for the 
ARP matrix. However, we still do not have a feeling of what the effect of Marp is when 
applied to an arbitrary Bloch vector. For getting it, we will calculate its eigensystem 
MARpcij = Ajaj. The eigenvalues and eigenvectors are 



Ai 



ai 



with 



1 





A, 



a2(A) 



/ cos((^(A)) \ 
sin((^(A)) 



v 







^(A) 



-X(A) + 



/ 

T 
2" 



A, 



«3(A) 



/ cos((^(A) + 
sin((p(A) + f ) 




(24) 



(25) 



We see from the set of eigenvalues that Marp is a rotation matrix of angle vr. The 
rotation axis is 03. This is the first relevant conclusion of our analysis: the effect of an 
ARP on a Bloch vector is equivalent to that of a tt rotation about an axis contained in 
the equatorial plane of the Bloch sphere. This by no way means that the Bloch vector 
actually performs this rotation. As we have discussed above, in the frame R' , the Bloch 
vector precesses around the control vector as this one goes from w' to —w'. To get the 
motion in the frame i?, we still need to compose that motion with the rotation about 
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w involved in the R to R' frame change. Indeed, the trajectory of the Bloch vector 
is much more complex than a vr rotation. However, if we take snapshots of the Bloch 
vector right before and right after an ARP, a vr rotation about an axis contained in the 
equatorial plane of the Bloch sphere links the two pictures. 

It is important to note that the rotation axis not only depends on parameters of 
the ARP pulse, such as 0(±T/2) or the time profile of VL{t) and It also depends 

on the transition frequency of the two-level system through x- We have highlighted this 
dependence by stating explicitly in ( l25l) that x ^^^id, hence, ip and are functions of 
A. This becomes relevant when dealing with an inhomogeneously broadened ensemble. 
In such a case, each frequency class experiences a vr rotation about a different axis, 
although all of these axes are contained in the equatorial plane of the Bloch sphere. For 
a given frequency class, the specific orientation of the rotation axis is determined by 
X(A). 

2.1.4- Comparison with a it pulse The ARP feature of being equivalent to a vr rotation 
about an axis contained in the equatorial plane reminds us of a vr pulse. In the latter 
case, if the external field is perfectly tuned to the two-level system, the Bloch vector 
performs a vr rotation about an axis contained in the equatorial plane. The differences 
between an ARP and a perfectly tuned vr pulse are, first, the rotation axes (except if 
the parameters for the ARP are especially chosen), and second, the fact that the Bloch 
vector does perform the rotation in the case of the vr pulse. 

Aside from that, the main dissimilarity between an ARP and a vr pulse rises when 
it comes to consider a spectral distribution of two-level systems. The dependence of 
these pulses on the detuning is quite different. In the case of an ARP, we can see from 

that the rotation axes for a range of A values span out on the equatorial plane. 
Nonetheless, the rotation angle is vr for all the two-level systems, despite their detuning. 
On the other hand, in the case of a vr pulse, the rotation axes for a range of A values fan 
out on a vertical plane, let us say, the plane u'w'. Moreover, the rotation angle is not vr 
for everybody. It is rather given by the expression Jinitf + A'^y/^dt, that reduces to 
[1 + (A/r2)^]^/^vr for a square pulse. The point is to analyse how one and the other kind 
of dependence on A impacts on the result of a sequence containing either an ARP or a 
vr pulse. 

Photon or spin echo: 

One of the most current uses of a vr pulse is the photon or spin echo, where the point 
is to rephase an inhomogeneous distribution of two-level systems. Starting from a set of 
Bloch vectors oriented along, let us say, the v' direction, the echo sequence, consisting 
in a waiting time r, a vr pulse and another waiting time r (symbolized r-vr-r), is ideally 
intended to yield — v'-oriented Bloch vectors. For the application of such a sequence to 
a detuned two-level system, we will define the rephasing error as the angle between 
the equatorial component of the final Bloch vector for the case A/f2 -C 1 and the final 
Bloch vector of the atom in resonance with the field. It can be shown that \e.,r\ scales 
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as (A/fi)^ (the exact expression for a square pulse is = sin(2Ar)(A/f2)^). 

Due to the similarities between the ARP and the tt pulse, it is valid to ask oneself 
whether the former can be used as a rephasing pulse in an echo experiment. Rephasing, 
if it occurs, will not necessarily happen with the Bloch vectors parallel to —v'. The 
direction of the rephased Bloch vectors will rather be the one the resonant vector takes 
at the end of the sequence. In any case, simple arguments immediately tell us that 
an ARP in an echo sequence does not seem a clever choice: As the orientation of the 
Bloch vector in the equatorial plane is a crucial issue in the rephasing experiment, an 
axial operation whose axis is variably oriented in the equatorial plane, as an ARP, will 
evidently not work. Let us now consider quantitative arguments. For the sequence 
r-ARP-r, we define £:arp in the same way as we did for a vr pulse. Simple calculations 
yield that £:arp scales as A^/r, where r is the mean value of the chirp rate (pit) (the 
equality £:arp = A^/r is obtained for a linearly chirped ARP). Due to the adiabatic 
condition (fT6|) . we deduce that I^arpI ^ i.^/^Y ~ I^ttI- We conclude that an ARP is 
unsuited for substituting a vr pulse in a regular echo sequence. At least, it will be much 
less efficient than the standard tt pulse. 

Population inversion: 

Another purpose for which vr pulses and ARPs are often used is population inversion. 
Let us compare the performances of both pulses. In this case, we will define e as the 
angle between the final Bloch vector and the desired vertical direction. For an ARP we 
have £arp = since, as we have seen, the ARP can be viewed as a tt rotation about 
an horizontal axis as long as fl22|) and f|T5l) are satisfied. As regards the tt pulse, for 
small A/n we get e.,^ ~ 2A/VL. It is clear, then, that the choice for population inversion 
purposes will be an ARP rather than a vr pulse. 

2.2. Rephasing by two adiahatic rapid passages 

2.2.1. The rephasing sequence matrix As we have seen in the previous section, an ARP 
is not able to rephase an inhomogeneously broadened ensemble of two-level systems 
unless Pinh is small enough or the chirp is fast enough to fulfil Tf^^^Jr ^1. It is 
known [Tl[2l[T8l|26], however, that two ARPs can produce the rephasing without any 
stringent condition other that the fulfuilment of (ITSl) and ( l22l) for the individual ARPs. 
In this section, we will analyse how the rephasing by two ARPs happens. 

Once again, we prefer to make use of the matrix representation of the process. 
The sequence to be considered is the following: the coherences are created at t = 0, the 
system evolves freely during a time ri, a first ARP (henceforth called ARP A) is applied 
between tx — Tx/2 = ti and tx + Tx/2, the system evolves freely again during a time r2, 
a second ARP (henceforth called ARP B) is applied between t-Q — T-q / 2 = ti + T2 + 
and t-B + 7b /2, the system evolves freely during a time rs. We read the coherence at 
instant t = ti + T2 + + Ta + Tb. The matrix associated to the sequence is 



L = F(r3)MARPBi^(r2)MARPAi^(ri), 



(26) 
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with 



F(r) 



v 







cos{uJabr) 








(27) 



the matrix for the free evolution time r in frame R, and 

Marpj = {Tj/2)r' [C2 {Tj/2)]-' Uj {Tj/2, -Tj/2) {-Tj/2) C^j {-Tj/2) , (28) 

the matrix for the ARP J = A, B. We have included de subindex J in Ci to indicate 
that for each ARP, the matrix in is to be calculated using the parameters Uq and 
of ARP J, henceforth noted Uq^j and (pj (subindexation is made extensive to all other 
ARP-dependent quantities). As for C2, the subindex is unnecessary since we consider 
that ([23D is valid for both ARPs. 

The matrix computation of L in (|26l) yields 

^ cos a — sin a ^ 
L = sin a cos a , (29) 



sm a 




cos a 




1/ 



with 



a = Uabin -T2 + n) - Xb(Ab) + Xa(Aa) + M-Tb/2) + MTb/2) 
-0a(-Ta/2)-0a(Ta/2). 

= Uabin -T2 + T3) + 2 [M^b) " ^a(Aa)] (30) 

We see that L is a counterclockwise rotation about w of an angle that depends on 
parameters of the ARPs but also on the atom transition frequency. For the sequence 
to succeed in rephasing an inhomogeneous distribution, we need to get rid of the 
dependence on Uab- This can be achieved by setting 

T3 = T-2 - n, (31) 

XA(Wa6 - C^O.a) = XB(Wa6 - Wo,b) ^ ^^ab ^ '^inh, (32) 

which leaves 

a = M-Tb/2) + MTb/2) - M-Ta/2) - MTa/2). (33) 

Fulfilling fl3T]) is trivial: it is enough to set T2 > ti and to impose the rephasing instant 
a time T2 — ti after the second ARP. As for ( l32l) . the simplest way to satisfy it is to 
use identical amplitude time profiles and identical frequency chirps for both ARPs. The 
phases need not be identical, though. Actually, it is the relative phase between ARPs 
both at the beginning and at the end of the pulses that determines the rotation angle 
of L or, in other words, the phase of the rephased coherence. If the phases of the ARP 
fields are also equal, which makes the second ARP just a time-shifted version of the 
first, L is the identity matrix. Another way of retrieving the identity matrix is to design 
the ARPs such that (f)j{Tj/2) + (j)j{—Tj/2) = 2n7T, with n an integer. 
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ARP A ARP B 




Figure 2. (Colour online) Field profile for two linearly chirped ARPs of constant 
n. The chirp rate, chirp range and field amplitude is the same for both ARPs. (a) 
Both ARPs have the same phase, i.e. (?!)a(±T/2) = (/)b(±T/2). This sequence gives 
phase preservation rephasing. (b) The ARPs have different phases, i.e. ^^(±T/2) ^ 
(/)b(±T/2). This sequences adds a phase factor to the rephased coherences. 

2.2.2. Phase preservation The above point about relative phases introduces a 
difference between rf and optical rephasing. Current rf technology permits accurate 
and arbitrary control not only of the amplitude but also of the phase of the rf field. 
The production of two identical pulses, identical both in their envelope and carrier, as 
depicted in figure |2]^a), is straightforward. Spin rephasing with such ARPs preserves 
the phase of the initial superposition state (L = I). On the other hand, the phase of 
an optical field is trickier to control. If the optical ARP pulses are to be produced from 
acousto-optical modulation of a monochromatic cw laser, the acousto-optic modulator 
allows one to control the relative phase between pulses, as long as the laser field stays 
coherent through the rephasing sequence. So the same situation as for the rf case is 
found. If the laser coherence time is shorter than the rephasing sequence, no phase 
control can be applied and the phase of the rephased coherence will be random. This is 
the situation sketched in figure [2t^b). 

For the processing or storage of quantum information encoded in the phase of the 
carrier photons (i.e. time-bin qubits), the above considerations become relevant. 

Having said that, we must prevent the reader from thinking that phase 
randomization can be avoided by preferring storage in the spin coherences. Acces 
to those coherences, either by conversion from optical ones or by EIT, involves the 
application of additional external light pulses. The random phase of those pulses results 
in phase randomization of the retrieved qubit as well. 

2.2.3. Geometric interpretation The rephasing process by two ARPs is easily 
understood by means of a simple geometric picture. We will go through the steps 
of the rephasing sequence analysing their action on the Bloch sphere. We will content 
ourselves with a 2D picture of what happens in the equatorial plane. In other words, we 
will only consider the effect of these operations on the projection of the Bloch vector onto 
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the uv plane. As we concluded in section 12. 1.3^ the effect of one ARP can be mapped in 
the Bloch sphere to a vr rotation about an axis contained in the equatorial plane oriented 
at a counterclockwise (positive) angle v^(A) = 1/2 [— x(A) + (— T/2) + (T/2)] from 
the axis i) (see (12^ and (!25|) ). The 2D version of this operation is a reflection in the 
mentioned axis. We will restrict to the case of two totally identical ARPs. Then, the 
rotation (3D) or reflection (2D) axis is the same for both ARPs. We will note this axis 
s, as shown in figure |3l 

Figure [3] shows, with very simple elements, how rephasing by two ARPs works. 
We will start from a Bloch vector Bq initially oriented along v (figure [3]^a)). From the 
geometrical point of view, the rephasing sequence consists in (see figure [3l): 

[i) The first free evolution period: A rotation of 9i = uJabTi- This leads from Bq to Bi 



(a) The first ARP: A reflection in axis s. This leads to B2 (figure [3](c)). 
(Hi) The second free evolution period: A rotation of 62 = uJabT2- This leads to B3 



{iv) The second ARP: A second reflection in axis s. This leads to B4 (figure [3]J^e)). 
(u) The third free evolution period: A rotation of 6*3 = 62 — 61. The Bloch vector regains 
its initial position Bq (figure |3]^f)). 

At instance [iv), we see that, to recover the initial vector Bq, all we need to do is to let 
B4 describe an angle 62 — 61. That is equivalent to condition fl3T]) . compulsory for the 
rephasing to take place, and is exactly what instance (v) is about. 

In figure [31 our choice of the angles ip, 61 and 62 has been totally arbitrary. Had we 
chosen a different set of angles, the result would have been the same. As s, 61 and 62 are 
determined by Uab, this proves that every Bloch vector in the inhomogeneous broadening 
regains its position of departure. This guarantees rephasing. In other words, we see from 
figure [3] that L\2d, the restriction of L to the equatorial plane, satisfies L\2d = I2x2- 

Regarding the Bloch vector vertical component, which we have neglected so far, 
the two 3D TT rotations combine to leave it unaffected. Hence, L = 1. 

If we consider ARPs with different phases (non- vanishing relative phase), the 
geometric interpretation becomes a little trickier because two different axes sa and 
sb are involved. Anyway, it is not hard to convince oneself that L is a rotation about 



w of angle 2{ipB - ^a) = 0b(-7'b/2) + 4>b{Tb/2) - 0a(-Ta/2) - 0a(Ta/2), just as in 



2.2.4- Comparison with ir-pulse rephasing As we stated in section 12.1.41 the error 
associated to rephasing by a vr pulse scales as (A/f2)^. The same behaviour is found 
for e2n, the error corresponding to a rephasing sequence with two vr pulses, either of the 
same or opposite rotation axis (the errors for these two cases differ at higher order of 



(figure [3]^b)). 



(figure [3t^d)). 



(EHD. 



A/n). In fact, 

£2. = [(sin(2Ari) - sin(2A(r2 - n)) - 2sin(Ar2) + 2sin(A(r2 - 2ri))] 




(34) 




Figure 3. (Colour online) (a-f) 2D geometrical representation of rephasing by two 
identical ARPs. The plane uv is the equator of the Bloch sphere. The initial vector Bq 

(a) describes a first rotation of angle 9i = iOabTi during the first free evolution period 

(b) . Then, the first ARP applies a reflection in axis s leading from Bi to B2 (c). A 
second rotation of angle 62 = uJab'T2 follows, corresponding to the second free evolution 
period (d). The Bloch vector gains position B3. From there, the second ARP takes 
it to B4 by means of a reflection in axis s (e). Finally, the last free evolution period 
applies a rotation of angle ^3 = L0ab{T2 — rl) that leaves the Bloch vector at its initial 
position (f). (g,h) Amplitude (g) and frequency (h) time dependences of the rephasing 
field. 

On the contrary, rephasing by two ARPs is A-independent as long as ( l22l) . ( |T5l) and 
f l32|) are fulfilled. Therefore, its associated error, e2ARP, is zero. 

For an experimental comparison of ARP- and 7r-pulse-based rephasing, see reference 

m- 

3. Experiments 

We have verified experimentally the rephasing time of the echo given by f l3T]) and the 
phase preservation of the initial superposition state ( 133|) with two identical ARPs. 

The material system used was the rare-earth ion-doped crystal Tm^+:YAG (dopant 
concentration: 0.1 at.%), of interest in quantum storage applications [SGj - Hl] . We have 
focused on rf ARP since accurate field control is available. A suitable rf spin transition 
is obtained in Tm^"*":YAG by ground level splitting under the application of an external 
magnetic field of a fraction of Tesla. This transition is inhomogeneously broadened due 
to the slightly different lattice environment seen by each Tm^"*" impurity. A convenient 
optical transition at 793 nm is used for initializing the system and probing its state. 
The technique used in the experiments was optically-detected NMR [iSl - HS] . 

The optical aspects of the setup have been described extensively in [37] and |10] . 
Basically, the light beam, emerging from an external cavity diode laser, is amplitude 
and phase-shaped by acousto-optic modulators, driven by a high sample-rate arbitrary 
wave form generator (AWG, Tektronix AWG5004). The crystal is cooled down to 1.7 K 
in a liquid helium cryostat. The static magnetic field is generated by superconducting 
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coils and oriented as in [37]. The rf ARP field is supplied by a 10-turn, 20-mm-long, 
10-mm-diameter coil oriented along the light pulse wave vector (the [110] axis of the 
cubic crystal lattice). The crystal sits at the coil centre. The rf signal, generated by 
the AWG, is fed to the coil through a 500- W amplifier (TOMCO BTOOSOO-AlphaSA) 
and a ~ 600-kHz-bandwidth resonant circuit. The amplitude of the static magnetic 
field (~ 0.5 T) is chosen such that the spin transition frequency matches the rf circuit 
resonance known to be close to 14 MHz. The inhomogeneous broadening of the spin 
transition in the sample was measured to be Finh ~ 500 kHz. 

The ideal initial state right before the application of the ARP-based rephasing 
sequence is that all the Bloch vectors of the atomic ensemble are aligned in the equatorial 
plane of the Bloch sphere (w', for example). We reach this situation in two steps, 
identical to those described in [1]. First, we make use of the 3-level A system of 
Tm^+:YAG under static magnetic field [36] to deplete one of the spin sublevels (noted 
level a) and fully populate the other one (noted level h) through optical pumping. This 
way we set B = wl. Then, the aim is to rotate that vertical Bloch vector to the 
equator. The usual way of achieving that is to apply a 7r/2 pulse. However, our limited rf 
power does not allow 7r/2 pulses short enough so that their bandwidth would cover Pinh- 
Therefore, we turn to an ARP pulse that is interrupted at half its way. The frequency is 
chirped from 0(— T/2) to 0(0). We will call such pulse an adiahatic half passage (AHP). 
If we neglect the inhomogeneous broadening and if the final frequency of the AHP (or 
the central frequency of the interrupted ARP) is tuned to the sublevel transition, the 
AHP (of positive chirp) turns the Bloch vectors from vol to n/, which is the desired 
situation. However, because of the inhomogeneous broadening, the Bloch vectors at the 
end of the AHP fan out typically from [VL{Q)u' - {Ti^^/2)w'] to [VL{{))u' + (rinh/2)ii;'] , 
as shown in figure Ht^a). This depart from the desired orientation ill will not be relevant, 
as we will see, since two identical ARPs recover the initial state of any Bloch vector 
(L = I). From this starting point, we execute the rephasing sequence described in 
section [2.2.11 At the end of the sequence, we apply a reversed AHP (rAHP), that is, an 
AHP where the frequency is chirped in the opposite sense between the same values of 
the previous AHP. If L satisfies L = I, the rAHP should bring all the Bloch vectors to 
wl, which was the state before the AHP. In fact, this is the main source of information 
in our experiments: the comparison between the states right before the AHP and right 
after the rAHP. 

To be more precise about the AHP and rAHP, their phase characteristics for a 
linear chirp are summarized as 



0AHp(t) = rtV2 and - T/2 < t < (35) 

0rAHp(^) = -rt^/2 and < t < T/2 (36) 

For probing the mean state of the system at any time of the experiment, we measure 
the intensity of a weak probe pulse tuned to the optical transition from state a to the 
upper state of the A system. We assume the probe pulse is weak enough not to alter the 
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Figure 4. (Colour online) (a) Span of Blocli vectors at the end of an AHP 
appied to an inliomogeneously broadened ensemble of atoms. The vector lying 
along it' corresponds to the atom in resonance with the last frequency of the AHP 
chirp (the central frequency of the interrupted ARP). The vectors in the upper 
([f](0)M' + (rinh/2)u>1/[l^(0)2 + (ri„h/2)2]V2) and lower ([l](0)M'-(ri„h/2)u>']/[l^(0)2 + 
(Finh/S)^]^/^) ends are the most detuned ones, at A = ztFinh/S. (b) Scheme of the 
plane of action of the rAHP control vector when an additional constant phase 5 is 
included in the field of the rAHP: (?l>rAHP = — + S. This plane is rotated an 
angle —S with respect to the plane u'w', where the initial AHP control vector acts 
((/)AHP = rt^/2). 

distribution of population between the levels. In such a case, the measured intensity / is 
a function of the population of level a as / = lQexp{—kpaa), with Jq the beam intensity 
before hitting the sample, k some positive constant determined by the opacity of the 
sample and paa the ensemble average of paa- Our conclusions will be drawn from the 
comparison between / before the AHP, /i, and / at after the rAHP, If. 

We performed two experiments. In the first, we aimed at testing the preservation 
of the initial superposition state, i.e. L = I when the two ARPs are identical. In the 
second, we verified the condition (13T]) : = T2 — Ti, which gives the rephasing time of 
the echo. 

3.1. State preservation: Validity of L = 1 for identical ARPs 

For the particular ARPs used in this experiment, the rf field amplitude at the centre of 
the circuit resonance yielded a Rabi frequency f2niax/(27r) = 141 kHz. The rf frequency 
was chirped linearly from lower to higher frequencies through a range of 4 MHz during 
100 fis. As a consequence of its resonance profile, the rf circuit modulates Q{t) even 
though the input rf current amplitude was kept constant throughout the ARP: Q{t) 
follows the circuit profile as the frequency is swept. This ensures fulfilment of the far 
off- resonance initial and final conditions f l22|) . To monitor the fulfilment of the adiabatic 
condition ( |T5|) . let us define the quantity C{t) as the left-hand side of (fT5l) . Using the 
parameters described above, ( is plotted as a function of t in figures M^a) and Mjo) 
for two different values of A, at the centre and edge of Pinh- We see that C < 1 is 
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(a) A = 

Q l{2n) = 141 kHz 



(b) A/(27t) = 250 kHz 
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Figure 5. (Colour online) Validity of the adiabatie approximation: ( is defined 
as the l.h.s. of plotted as a function of t for the experimental parameters 

(/>{±T/2)/(27r) = ±2 MHz, T = 100 /iS, (I)/{2tt) = 0.02(t//zs)2. Condition (Hi]) is 
not satisfied at every instant of the ARP for frequencies at the edge of Finh- (a) 
A = 0, f7max/(27r) = 141 kHz: C < 1 throughout the ARP. (b) A/(27r) = 250 kHz, 
f^max/(27r) = 141 kHz: ( can take values higher than unity. However, C > 1 during 
just 2% of the ARP. (c) A = 0, r2max/(27r) = 262 kHz. (d) A/(27r) = 250 kHz, 
f^max/(27r) — 262 kHz. The behaviours in (c) and (d) are similar to those in (a) and 
(b), respectively, except that ( takes lower values. 



satisfied by the experimental parameters at the centre of Finh (figure [5]^a)). However, 
in the case A/(27r) = 250 kHz (figure [5](b)), C < 1 is violated during a short period 
that amounts to just 2% of the ARP. In any case, we will see that this transient failure 
of the adiabatie character of the passages for some atomic frequency classes does not 
invalidate the conclusions drawn from the experiment. 

To ensure the preservation of the initial superposition state, i.e. L = I, quantum 
process tomography [19] would be the proper experiment to carry out. Here, we will 
content ourselves with simpler approximate tests. 

The first and simplest test is to check if Jf = Ji. Figure [6](a) shows the intensity 
profile of the probe beam as the rephasing sequence takes place. We will focus on the 
curve labelled "0". At t = —50 /xs, the intensity Ii corresponds to the transparency 
value (level a is depleted by previous optical pumping). The intensity decreases slowly 
at a first stage, but close to t = the intensity fall speeds up. It coincides to the short 
period where f3 is far off the vertical axis. At t = 0, / attains a value corresponding to 
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t{\xs) 

Figure 6. (Colour online) (a) Probe intensity profile through a double ARP rephasing 
sequence. An AHP (rAHP) is applied before (after) the sequence. The parameters for 
both ARPs are fJniax/(27r) = 141 kHz and (j){t)/{2'K) = 0.02(t/^s)2. Spin T2 in this 
sample has been measured to be 510 (xs. Several phases of the rAHP are tested. 
For each curve, the rAHP phase is set to — (/)arp + 5. The labels indicate 5/{2it). (b) 
Experimental (symbols) and theoretical (line) ratios between final and initial intensities 
as a function of the phase of the rAHP. 

pbb = paa- After a free evolution time of 10 /is (the probe beam is turned off during this 
and the remaining free evolution intervals), the first ARP starts. The intensity varies 
slowly during the initial and final regions of the ARP, but it displays a peak at half its 
way. The origin of this peak is the fact that the Bloch vectors fan out in a circle section 
contained in the u'w' plane at the end of the AHP, as depicted in figure Ht^a). If all 
Bloch vectors were aligned at t = 0, no peak would be observed. Another point is that 
the intensity should be the same on both sides of the peak (see figure [7] as an example). 
We estimate the intensity difference is due to imperfect matching of the rf frequency to 
the centre of the inhomogeneous broadening. After a second free evolution time of 20 
/is, the second ARP is apphed. In this case, we observe a dip instead of the peak, of 
the same origin than that of the latter. The last free evolution interval of 10 /is takes 
place. Finally, the rAHP is applied. The intensity increases to a value /f close to the I\. 

The experiment shows that If < However, we cannot yet conclude that the 
reason for Jf lower than expected is imperfect rephasing (L 7^ I). As a matter of fact, 
the observed If is totally compatible with the decrease induced by finite spin coherence 
lifetime T2. The latter was previously measured to be ~ 510 /is (results not shown). 
Therefore, as far as this first test is concerned, rephasing is well-achieved by the double 
ARP sequence. The transient failure of the adiabatic condition for some frequency 
classes described above does not seem to harm the rephasing efficiency. 

The second test consists in monitoring If as the phase of the rAHP is varied from 
its optimal value given by (|36|1 . If a phase — + 5 is assigned to the linearly chirped 
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rAHP, the control vector during the rAHP will no longer act in the vertical plane u'w' . 
It will act in another vertical plane rotated about w' an angle —5 from the original one 
(see figure 111(b)). As a consequence, this phase-shifted rAHP will be most suited for 
driving a Bloch vector that is contained in the same plane at the end of the rephasing 
sequence. In other words, a rAHP of phase —rt^/2 + 5 will best fulfil its aim of returning 
the Bloch vectors to the vertical axis if L is a rotation matrix as in fl29|) with a = —6. 
Thus, this experiment will allow us to check if a is indeed zero in the case of a rephasing 
sequence with identical ARPs. 

The curves obtained with different values of S are exhibited in figure [6]^a). The 
label of each curve corresponds to 6/{2tt). We observe that If varies as a function of 6, 
displaying highest values around 6 = and lowest around 6 = n. This means that when 
5^0 (5 ~ vr), the rAHP guides the Bloch vectors mainly to the w' {—w') direction. 
The comparison of the experimental ratio /f//i with the predictions based on the theory 
in section [2.2.11 is presented in figure [6]^b). We observe that the maximum is somewhat 
shifted to the left of the expected value 6 = 0. The minimum ratio is also not as low as 
predicted by the theory. These differences might be due to the slight mismatch between 
the centre of the frequency chirp range and the centre of Pinh and/or to the transient 
failure of the adiabatic condition for some frequency classes. In any case, it is clear that 
the overall behaviour is compatible with a sequence matrix close to the one obtained 
for L with a = (0.05 ± 0.1)27r. 

3.2. Rephasing time condition: Validity of = T2 — ti 

For this experiment, the rf circuit was improved to reach a higher Rabi frequency: 
^max/(2vr) = 262 kHz. This allowed us to better satisfy the adiabatic condition ( ITSjl . 
as shown in figures [TST c) and ITST d). The remaining parameters of the ARPs were kept 
identical to those in the previous experiment. 

We studied the effect on /f of varying ti while keeping T2 constant. For each 
rephasing sequence, we adapted rs to fulfil = T2 — ti (equation ( 13T|) ). The results 
are presented in figure [71 Both panels contain the same data (just the way the data 
are displayed changes). First of all, we observe that the quality of the curves is better 
than in figure |6l(a): the intensity level on both sides of the peak (or dip) are better 
balanced, which is a signature of better matching of the rf frequency to the centre of 
the inhomogeneous broadening. Second, the curves look very similar to one another, 
except, of course, for the positions of the ARPs. The resemblance is better appreciated 
in figure 17(b), where we clearly see that /f is the same for all curves. Here again, the 
value of If (< /i) is compatible with the spin decoherence and with the value of T2. 
Therefore, we conclude that, as far as condition fl3T|) is fulfilled, the performance of the 
rephasing sequence is optimum, independently of the particular values of ri and T2. 
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Figure 7. (Colour online) ARP rephasing experiments with varying ti. T3 is adjusted 
to satisfy T3 = r2 — ti. (a) Stacked curves, (b) Superimposed curves. 



4. Summary 

We have thoroughly analysed how rephasing of optical or spin coherences by two ARPs 
in a double echo scheme works. In the frame of the Bloch sphere formalism, we have 
theoretically developed a matrix treatment for the sequence, where the resulting matrix 
L is the product of the matrices associated to each building block. We first derived the 
matrix for a single ARP, which turned out to be a rotation matrix of an angle vr about an 
axis contained in the equator of the Bloch sphere. The particular orientation of this axis 
within the equator depends both on specific parameters of the ARP and on the atomic 
frequency. The latter dependence is the reason why one single ARP cannot manage to 
rephase an inhomogeneously broadened distribution of atomic coherences. A sequence 
involving two ARPs, however, is able to achieve the rephasing. In our matrix approach, 
this is understood as L = I, as derived in our calculations as long as some particular 
conditions are satisfied. We identify these key conditions. One can be satisfied very 
easily: it is enough to wait a longer time between the two ARPs (time interval T2) than 
between t = and the first ARP (time interval ri). Then, the rephasing takes place 
a time T2 — ti after the end of the second ARP. The second condition implies that the 
Rabi frequencies and frequency chirps of both ARPs must be identical. The rephasing 
process can be easily explained with the help of a very simple geometrical interpretation 
that is drawn from our matrix treatment. 

We have also analysed the capability of the sequence to preserve the initial state 
phase, assuming the two conditions just mentioned are satisfied. This is of particular 
importance for quantum memory applications. We have found that phase preservation 
is assured if the optical or rf fields of both ARPs have exactly the same time-varying 
phase, meaning that the second ARP field must be a time-shifted copy of the first. This 
opens a discrepancy between optical and rf rephasing, since the ability of controlling 
the field phase is different in the rf and optical technologies. In addition, we have 
provided a quantitative comparison between rephasing by vr pulses, as in standard echo 
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experiments, or by ARPs. The rephasing efficiency is superior in the latter case. 

We have verified experimentally the two rephasing conditions in the rare-earth ion- 
doped crystal Tm^+:YAG. Optically detected NMR experiments evidenced the rephasing 
capabilities of rf ARP rephasing sequences. If the finite lifetime of the coherence (T2) 
is taken into account, the experiments compare rather satisfactory to the theoretical 
predictions. 
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